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Abstract. We present a new method to obtain spatio-temporal information from aggregated data of sta- 
tionary traffic detectors, the "adaptive smoothing method". In essential, a nonlinear spatio-temporal lowpass 
filter is applied to the input detector data. This filter exploits the fact that, in congested traffic, perturba- 
tions travel upstream at a constant speed, while in free traffic, information propagates downstream. As a 
result, one obtains velocity, flow, or other traffic variables as smooth functions of space and time. Appli- 
cations include traffic-state visualization, reconstruction of traffic situations from incomplete information, 
fast identification of traffic breakdowns (e.g., in incident detection), and experimental verification of traffic 
models. 

We apply the adaptive smoothing method to observed congestion patterns on several German freeways. 
It manages to make sense out of data where conventional visualization techniques fail. By ignoring up to 
65% of the detectors and applying the method to the reduced data set, we show that the results are robust. 
The method works well if the distances between neighbouring detector cross sections do not exceed 3 km. 

1 Introduction 

During the last decades, traffic dynamics and pedestrian flows have been intensively 
studied. Regarding the observed phenomena and simulation approaches we refer 
the reader to some recent reviews []IJ£||3|,|J. Presently, scientists are more and more 
getting interested in detailed empirical studies. Reasons for this are the availability 
of better data and the need to verify and calibrate models. Both, theoretical 
and empirical p|,|9|, p!0| , pT| . p^| . p!3| , p!4| . p3| . pTj| , p!7| studies indicate that the phenomenology 
of congested traffic is more complex than originally expected. There seems to be a 
rich spectrum of traffic states [p|,|6|, [T2| , [T3[1 , hysteretic or continuous temporal or spatial 
transitions among them [j7j pl] , [i^ . [i^|i^[l9|1 , and fluctuations or a erratically appearing 



dynamics play a significant role p3|,p!6|j20|,|2T|,|22| . In order to make sense out of this, it 



is increasingly important to have suitable ways of data processing, in order to extract 
the information relevant for scientific investigations or specific applications. Here, 
we will propose a three-dimensional data-evaluation method allowing to visualize 
the spatio-temporal dynamics of traffic patterns along freeways. 

The developed "adaptive smoothing method" filters out small-scale fluctuations 
and adaptively takes into account the main propagation direction of the information 
flow (i.e. the dominating characteristic line), which have been determined by means 
of a spatio-temporal correlation analysis in other studies p3| . The temporal length 



scale of the smoothing procedure can be as small as the sampling interval, while the 
spatial length scale is related to the distance between successive detectors, which 
can be up to 3 kilometers long. 

By "filter" we just mean a transformation of the data with specific properties. 
Here, we use a spatio-temporal lowpass filter, i.e., only (Fourier) components of 
low frequency can pass the filter, while high-frequency contributions are considered 
as fluctuations and smoothed out. One particular feature of our filter is that it is 



nonlinear and adaptive to the traffic situation in distinguishing free and congested 
traffic, as the propagation direction of perturbations differs. The results are three- 
dimensional visualizations of traffic patterns, which are quite robust with respect to 
variations of the filter parameters and very helpful in obtaining a clear picture of 
the systematic spatio-temporal dynamics. 

Therefore, this method is suitable for the reconstruction of traffic situations from 
incomplete information, fast identification of traffic breakdowns (incident detection), 
and experimental verification of traffic models. First results support the phase dia- 
gram of traffic states occuring at bottlenecks |I7|J2^| , which, apart from free traffic, 
predicts pinned or moving localized clusters, spatially extended patterns such as 
triggered stop-and-go waves, and oscillating or homogeneous congested traffic, or a 
spatial coexistence of some of these states P,|7|,[l]]. The respectively occuring spatio- 
temporal pattern depends on the specific freeway flow and bottleneck strength, but 
also on the level of fluctuations, as these can trigger transitions from, for example, 
free traffic to localized cluster states PJT7|JT|1. We do not see sufficient support for one 



"generalized pattern" p5| that would always be observed when traffic flow breaks 
down. 



2 Description of the Method 

The adaptive smoothing method is a data processing method for obtaining traffic 
variables as smooth functions of space and time out of stationary traffic data. It has 
following heuristically motivated properties: 

1. In case of free traffic, perturbations (of, e.g., velocity or flow) move essentially 
into the direction of traffic flow More specifically, they propagate with a 
characteristic velocity Cf ree at about 80% of the desired velocity Vq on empty roads 
23| . Therefore, at locations with free traffic, perturbations with propagation 



velocities near Cf ree should pass the filter. 

In case of congested traffic, perturbations propagate against the direction of 
traffic flow with a characteristic and remarkably constant velocity c con g ~ — 15 
km JT2]]. With modern data analysis techniques, it has been shown that such 



propagation patterns persist even in "synchronized" congested traffic flow, where 
they are hardly visible in the time series due to a wide scattering of the data 
in this state |26|]. So, for high traffic densities or low velocities, the filter should 
transmit spatio-temporal perturbations propagating with velocities near c con g 
more or less unchanged. 
3. The filter should smooth out all high-frequency fluctuations in t on a time scale 
smaller than r and spatial fluctuations in x on a length scale smaller than a. 
The parameters r and a of the smoothing method can be freely chosen in a wide 
range (cf. Table [1]). 

Let us assume that aggregated detector data zf- are available from n cross sec- 
tions i at positions Xi where j G { j m in, • • • , jmax} denotes the index of the aggregation 
intervals. 

Usually, the aggregation interval 



At = t j -t j - 1 (1) 



is fixed (between 20 s and 5 min, depending on the measurement device). A typical 
value for German highways is At — 1 min. 

The components of z™- represent the desired aggregated quantities as obtained 
from cross section % during the time interval j. Typical examples include the average 
velocity Vij, the vehicle flow Qij, the occupancy O^, or some derived quantity such 
as the traffic density py. The input data z 1 ? can represent either averages over all 
lanes or quantities on a given lane. 

The adaptive smoothing method provides estimates z(x,t) for all locations x G 
[xi,x n ] between the positions of the first and the last detector, and for all times 
t G [t m - m , tmax]. Extrapolations (in the sense of a short-term traffic forecast) will also 
be discussed. Without loss of generality we assume x\ < X2 < ■ ■ ■ < x n , and a traffic 
flow in positive x-direction. 

The core of our "adaptive smoothing method" is a nonlinear filter transforming the 
discrete input detector data z 1 *- into the smooth spatio-temporal functions z(x ) t). 
To satisfy the first two requirements mentioned above, we write the filter as 

z(x, t) = w(z 

cong) ■2'free)'2'cong('£) ^) [l ^(-2-cong; % free)] ^ free 

t). (2) 

This is a superposition of two linear anisotropic lowpass filters z cong (x, t) and 
Zf ree (x, t) with an adaptive weight factor < w < 1 which itself depends nonlinearly 
on the output of the linear filters as discussed later on. 
The filter z cong (x, t) for congested traffic is given by 

71 Jmax 

z cong( x ,t) — "77 / T\ ^1 (pcongi^i ~ X,tj — tyz^, (3) 

Jv cong^ ? IJ ,_ , ._. . 

L — ± J — jmin 

with 



tj = t min + 3 At, (4) 

and the normalization factor 



Tl Jmax 



N cong (x, t) = ^ XI <^ong(Zi - X, tj - t). (5) 

* — 1 J — Jmin 

This normalization guarantees that a constant input vector z™ = Zq for all i,j is 
transformed into a constant output function with the same components: z(x, t) = zq. 
Notice that the normalization depends on both, location x and time t. 

The analogous expression for the lowpass filter Zf ree (x,t) for free traffic is 

Tl Jmax 

Z fTCC (x, t) = — — - ^ Yl ^ee(^ - X, tj - t)zfj (6) 

* — L J — Jmin 

with 

Tl Jmax 

N ilee (x, t) = Y Y ^(^i - X, tj - t). (7) 

* 1 J Jmin 
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Fig. 1. Visualization of the effects of linear homogeneous filters with the kernels 
<f)f Tee (x, t) and 0con g (£, t), respectively. The shaded areas denote the regions consid- 
ered in the calculation of a data point at (x,t). Triangles denote the mainly con- 
tributing input data sampled in free traffic, squares the ones sampled in congested 
traffic. 

The kernels (p CO ng(x,t) and (f){ Tee {x, t) of the linear homogeneous filters do the 
required smoothing and are particularly transmissible for perturbations propagat- 
ing with the typical velocities c con g and Cf ree observed in congested and free traffic, 
respectively (cf. Fig. [I]). 

If the propagation velocity of the perturbations was zero, the filter should only 
perform the smoothing. We chose the (non-normalized) exponential function 

<f) (x, t) = exp (-^j- ~ ■ (8) 

Instead, one could apply other localized functions such as a two-dimensional Gaus- 
sian. However, it turned out that the exponential had more favourable properties in 
our application. 

The linear filter for nonzero propagation velocities can be mapped to 4>o{x, t) by 
the coordinate transformations (cf. Fig. |TJ) 

x — x , t = t cong H = t free H . (9) 

Ccong Cf ree 

Thus, we obtain for the kernels of the linear anisotropic filters 

0con g (a;, t) = <f) (x', t' cong ) = <f) lx,t — J , (10) 

free (x, t) = (j) (x', 4 ee ) = 00 ( X, t ) . (11) 



Figure || shows the action of the filters z cong (x,t) and Zf Tee (x, t) for the velocity 
fields V COIig (x,t) and V iree (x, t). 




Fig. 2. Typical velocity fields Vf vee (x,t) (top) and V cong (x,t) (bottom) obtained by 
application of the filters Zf ree (x, t) and z cong (x,t) to traffic data of a section of the 
German freeway A9 South. 

Finally, we define the nonlinear adaptive weight function w(z cong , z ircc ) 6 [0, 1]. 
Obviously, we must have w ~ I for congested traffic, and w ~ for free traffic, so 
we need some a priori estimate of the traffic situation at the point (x,t). 

Congested traffic is characterized by a high traffic density and low average ve- 
locity. Since, in contrast to the density, the velocity can be directly measured with 
stationary detectors, we chose the velocity to determine the a priori estimate. Dif- 
ferent possibilities to estimate the velocity at point (x, t) are: 

— The measured velocity of the detector cross section whose position is 
nearest to x in the time interval j containing the actual time t, 

— the velocity V cong (x,t) as calculated with the "congested-traffic" filter z cong ac- 
cording to Eq. (0) with kernel (|T0| ) (assuming a priori congested traffic), 

— the velocity Vf ree (x, t) as calculated with the "free-traffic" filter Zf ree according to 
Eq. (|) with the kernel (0), 

— or some combination of the above estimates. 



The first way to estimate the velocity is subject to errors, if the typical length scale 
A of occuring stop-and-go structures is not larger than the distance Axi between 
two neighbouring detectors. At this point, it is crucial that propagating structures in 
congested traffic, especially stop-and-go waves, are very persistent. It has been shown 
]27| that they can propagate through freeway intersections or other inhomogeneities 



nearly unchanged, passing all perturbations of free traffic on their way (see, e.g., Fig. 
D). Therefore, whenever at least one of the estimates indicates congested traffic, the 
weight function should favour the filter for congested traffic. Specifically, we assume 



cong; -2-frcc 



w(V coag (x, t), Vf ree (x, t)) 



tanh 



V c -V*(x,t) 
AV 



where 



V*(x, t) = min V cong (x, t), V {ree (x, t) 



, (12) 



(13) 



V c and AV are parameters that can be varied in a wide range. 

If not explicitely stated otherwise, for all simulations in the following section 
we will use the parameters specified in Table |l| We will also show that the four 
parameters c cong , Cf ree , V c and Av can be varied in a wide range without great 
differences in the output. In this way, we show that the proposed adaptive smoothing 
method does not need to be calibrated to the respective freeway. One can take the 
values from Table [1] as a global setting. 

The smoothing parameters a and r have the same meaning and the same effect 
as in standard smoothing methods, e.g., Eqs. (22) and (23) in Ref. JT7 . 



Parameter 


Typical Value 


Meaning 


a 


0.6 km 


Range of spatial smoothing in x 


T 


1.1 min 


Range of temporal smoothing in t 


Cfree 


80 km/h 


Propagation velocity of perturbations in free traf- 
fic 


Ccong 


-15 km/h 


Propagation velocity of perturbations in con- 
gested traffic 


Vc 


60 km/h 


Crossover from free to congested traffic 


AV 


20 km/h 


Width of the transition region 



Table 1. Parameters of the "adaptive smoothing method" defined by Eqs. 
their interpretation, and their typical values. 



-(H) 



In summary, the proposed "adaptive smoothing method" is given by Eq. (|2|) with 
the nonlinear weight function (|T^), the filter (^) with normalization (|5|) and kernel 
fllpp for congested traffic, the filter (^|) with normalization (^) and kernel (|TT]) for free 
traffic, and the smoothing filter (|3|). Table |T| gives an overview of the six parameters 



involved and typical values for them. The adaptive smoothing method includes the 
following special cases: 

— Isotropic smoothing resulting in the limits Cf ree — > oo and c cong — ► oo (for practical 
purposes, one may chose Cf ree = c cong = 10 6 km/h); 

— only filtering for structures of congested traffic in the limit V c — > oo (for practical 
purposes, one may set V c ^> Vo with Vo being the desired velocity); 

— only filtering for structures of free traffic for V c = AV = 0; 

— consideration of the data of the nearest detector only, if a = 0; 

— application of the actual sampling interval of the detectors, if r = 0, Cf ree — > oo 
and c cong — > oo. (To avoid divisions by zero, zero values of a, r, and AV are 
replaced by very small positive values in the software). 

3 Application to German Freeways 

We will now discuss data from the German freeways A8-East and A9-South near 
Munich, and of the freeway A5- North near Frankfurt. In all cases, the traffic data 
were obtained from several sets of double-induction-loop detectors recording, sepa- 
rately for each lane, the passage times and velocities of all vehicles. Only aggregate 
information was stored with an aggregation interval of At — 1 min. We will use the 
following input data z 1 ^: 

— The lane-averaged vehicle flow 

i=i 

where n\ - is the vehicle count at cross section i during time interval j on lane I. 
The considered sections of the freeways have L = 3 lanes in most cases. 

— The lane-averaged mean velocity 

f« = E%^, US) 
i=i 

where is the average velocity at cross section i during time interval j on lane 
I. 

— The traffic density determined via the formula 

Qij 

pij = W ' 

as occupancies were not available for all freeways. 

Figure ^ shows an example of a complex traffic breakdown that occurred on 
the freeway A8 East from Munich to Salzburg during the evening rush hour on 
November 2, 1998. Two different kinds of bottlenecks were relevant, (i) a relatively 
steep uphill gradient from x = 38 km to x = 40 km (the "Irschenberg"), and (ii) an 
incident leading to the closing of one of the three lanes between the cross sections 
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Fig. 3. Complex congested traffic pattern on the German freeway A8-East from 
Munich to Salzburg during the evening rush hour on November 2, 1998. (a) Sketch 
of the freeway, (b) Plot of the spatio-temporal density p(x, t) using conventional 
smoothing (resulting for the setting Cf ree = c con g = 10 6 km/h). (c) Plot of the same 
data as calculated with the adaptive smoothing method. 



D23 and D24 between t = 17:40 h and t = 18:10 h. For further details, see Ref. [[L7 
Figure ||(c) shows the traffic density using the proposed data processing method 
with the parameters specified in Table |]. For comparison, (b) shows the result of 
conventional isotropic smoothing ignoring the speed of information propagation, cf. 
Ref. JT7| . While the structures in the free-traffic regions are nearly the same in both 
plots, the new smoothing method gives a better reconstruction of the congested 
traffic patterns. In contrast to the conventional method, our proposed method can 
resolve oscillations with wavelengths A comparable to or smaller than the distance 
Axi between two neighbouring cross sections as illustrated for the three stop-and-go 



waves propagating through the whole displayed section (A ~ 2 km). Notice that the 
conventional method generates artefacts for this example. 

The congestion patterns caused by the uphill gradient (the "wings" in Fig. |3](b), 
(c)) have larger characteristic wavelengths and are resolved, at least partially, also 
with the conventional method. 




Fig. 4. Time series of the empirical density, approximated by means of Eq. (|16|), for 
several cross sections. 



For comparison, we plot time series of some cross sections in Figure Let us 
consider more closely the jammed traffic patterns caused by the temporary bottle- 
neck associated with the incident. While the fluctuations are irregular at detector 
D22 (about 1.2 km upstream of the aforementioned incident), they grow to three 
stop-and-go waves further upstream, cf. the arrows in Fig. |5[ However, due to mea- 
surement errors and other errors in determining the density , the time series of the 
density is very noisy, especially that of cross section D20. The adaptive smoothing 
method suppresses this noise without suppressing the structure of the pattern. 

To check, whether there are significant patterns in the data of detector D20 at 
all, we plot the velocity data at D20 as well (Fig. |5|). In contrast to the density, 




Fig. 5. Time series of the empirical velocity data determined via Eq. (15) for cross 
section D20. 



the velocity is measured directly resulting in less noise. One clearly sees three dips 
corresponding to three stop-and-go waves which confirms the results of the "adaptive 
smoothing method". 




Fig. 6. Spatio-temporal density p(x, t) as in Figure |3|, but using only the data of 
cross sections D16, D18, D20, and D22 as input (cf. Figure f|). 



As a further check that the structures shown in Fig. |3](c) are not artefacts of the 
data processing, we tested the method for a reduced data set. Figure |j] shows the 
result using only the even-numbered detectors as input, whose time series are plotted 
in Fig. |j. The main structures remain nearly unchanged, particularly, the three stop- 
and-go waves propagating through the whole section. Notice that the whole plot is 
based on data of only four cross sections, and that the data are extrapolated in the 
upstream direction over a distance of nearly 2 km. 

We now apply the processing method to data from a section of the freeway A9- 
South near Munich. There are two major intersections II and 12 with other freeways. 
Virtually every weekday, traffic broke down to oscillatory congested traffic upstream 
of each of these intersections. (For details, see Ref. JT7|] ). 
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Fig. 7. Two regions of stop-and-go waves on the German freeway A9-South near 
Munich occurring upstream of the freeway intersections II and 12. (a) Sketch of the 
freeway, (b) Plot of the inverse l/V(x, t) of the spatio-temporal velocity V(x,t), us- 
ing conventional smoothing, (c) Plot of the same data using the adaptive smoothing 
method. 



Figure [7] (b) shows the inverse of the velocity, 1/V(x,t) for a typical congested 
situation. When plotting 1/V(x,t), the structures of congested traffic comes out 
more clearly than for the density. The method resolves small density clusters in the 
region 508 km < x < 510 km between about 8:30 am and 9:30 am which disappear 
around x = 508 km. These structures cannot be identified with the conventional 
smoothing method (see Fig. 0(c)), which just shows a hilly pattern. 

Figure 0(b) is an example for the spatial coexistence of homogeneous congested 
traffic (a short stretch around x = 510 km), oscillating congested traffic (around 
x = 509 km), and stop-and-go waves (around x = 507 km). This may be a three- 



dimensional illustration of the so-called "pinch effect" [|TT|,|Tj], but we do not observe 
a merging of narrow clusters to form wide moving jams. The narrow structures 
of short wavelengths rather disappear. This may be an artefact of our smoothing 
method, so that video data are required to get a more detailed picture. Simulations 
of this spatial coexistence, however, indicate that narrow clusters rather disappear 
than merge, namely when they do not exceed the critical amplitude in an area of 
metastable traffic 0. 

Notice that the bottleneck causing this congestion is located at about x = 510 
km implying that it is caused by weaving traffic and slowing down at an off-ramp, 
not by an on-ramp as is most often the case on German highways. 




Fig. 8. Same data and same processing method as in Figure 0(c), but with (a) a 
propagation velocity of c con g = —12 km/h instead of -15 km/h (b) c con g = —20 km/h, 
(c) with the crossover parameter V c = 40 km/h instead of 60 km/h. 



Now we demonstrate that the method is robust with respect to reasonable pa- 
rameter changes. Figure |8|(a) shows the result for an assumed propagation velocity 
of c cong = —12 km/h instead of -15 km/h. In plot (b), we assumed c cong = —20 km/h 
instead. In both cases, the results are systematically better than with the isotropic 
procedure. Since the propagation velocity of perturbations in congested traffic is 
always about c cong = —15 km ]!r2"| , p3| , pgpm , one can take this value as a global 
setting. 

Due to the comparatively high magnitude of the propagation velocity Cf ree , the 
related time shifts in the transformation for free traffic are small. Therefore, the 
method is insensitive to changes with respect to this parameter. It turned out that 
taking the isotropic limit Cf ree — > oo nearly gives the same results. 

The crossover parameters V c and AV can be varied within reasonable limits as 
well. As an example, Fig. |8|(c) shows the result after changing the crossover velocity 
from V c = 60 km/h to V c = 40 km/h. 

Finally, we applied the method to isolated moving localized clusters and pinned 
localized clusters observed on the A5-North near Frankfurt. This freeway is particu- 
larly well equipped with detectors, so that the typical length scales of the observed 
structures are always larger than the distance between neighbouring detectors. How- 
ever, we obtained nearly the same results when using only the even-numbered de- 
tectors or only the odd-numbered detectors as input. Figure |9] shows that even a 
further reduction of the information to only 35% of the detectors yields good results. 

4 Summary and Outlook 

We have proposed a new "adaptive smoothing method" for the three-dimensional 
visualization of spatio-temporal traffic patterns, which takes into account the char- 
acteristic propagation velocities observed in free and congested traffic. The method 
is robust with respect to variations of its parameters, so that it can be applied 
to new freeway sections without calibration. In principle, it would be possible to 
determine the parameters (such as the propagation velocity of perturbations) lo- 
cally (e.g., by means of a correlation analysis). However, the results would look less 
smooth and regular, as the small number of data to determine the local parameters 
would be associated with considerable statistical errors. Consequently, a large part 
of the variations in the local parameters would not reflect systematic variations of 
the parameters. So, both the local and global parameter calibration may produce 
artefacts, but it is advantageous to use global parameter settings. Moreover, there 
is empirical support for surprisingly constant propagation velocities Cf ree and c con g of 
perturbations in free and congested traffic. The parameters V c and AV are related 
to the transition from free to congested traffic and, therefore, can also be well de- 
termined. The spatial and temporal smoothing parameters a and r can be specified 
according to the respective requirements. Suitable parameters allow a good repre- 
sentation of traffic patterns even when the distances between successive detectors 
are about 3 kilometers. 

We point out that the suggested "adaptive smoothing method" itself can be ap- 
plied to calibrate the characteristic propagation velocities Cf ree and c cong . For this 
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Fig. 9. Freeway A5-North from Frankfurt to Kassel during the evening rush hour 
on August 7, 1998. (a) Sketch of the freeway, (b) Data reconstruction of the velocity 
using all detectors D10-D26. (c) Reconstruction using a reduced data set of only 6 
of all 17 detectors used in (b). 

purpose, the ranges a and r of spatial and temporal smoothing are chosen small. 
The optimal propagation velocities minimize the offsets in the propagation patterns. 

The aim of the "adaptive smoothing method" is to reconstruct the spatio- 
temporal traffic data from incomplete information as good as possible to allow a 
better understanding of the complex traffic dynamics. Potential applications are, 
for example, traffic state visualization, incident detection, or the experimental ver- 
ification of traffic models. Our method could be further improved by taking into 
account information about the traffic dynamics such as the continuity equation or 
a suitable equation for the average vehicle velocity as a function of space and time. 
It could, then, be used to determine short-term traffic forecasts. 
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